On the presence of mid-gap states in CaX^Og 
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Using exact diagonalizations of finite clusters with up to 32 sites, we study the Ji — J2 model on 
the 1/5 depleted square lattice. Spin-spin correlation functions are consistent with plaquette order 
in the spin gap phase which exists for intermediate values of J2/J1. Besides, we show that singlet 
states will be present in the singlet-triplet gap if J2/J1 is not too small (J2/J1 > 0.47). We argue 
that this property should play a central role in determining the exchange integrals in CaV-iOg 
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The interest in 2D frustrated magnets has recently in- 
creased with the discovery of the first 2D spin 1/2 system 
exhibiting a spin gap, namely CaV^g This material 
consists of V4O9 planes separated by calcium layers |Q, 
each vanadium atom being at the centre of a VO5 square- 
pyramid of oxygens. Vanadium having an oxidation num- 
ber 4+ in this compound, there is a single electron in the 
d-shell of each vanadium atom. These electrons behave 
as localized spins 1/2 coupled by some exchange inter- 
actions. The relative magnitudes of the various possi- 
ble exchange mechanisms cannot be easily deduced from 
quantum chemistry, and to get a better understanding of 
the magnetic properties of this compound, a direct com- 
parison of physical properties of model Hamiltonians to 
experimental results is necessary. 

The minimal model to describe this compound, namely 
the 1/5 depleted Heisenberg model with nearest neigh- 
bour interactions J\, has been studied in great details. 
Contrary to early results M, it is now clear that this 
model has Neel long range order and no spin gap. This 
was first proposed on the basis of exact diagonalizations 
and Schwinger boson mean field calculations Q and con- 
firmed by Monte Carlo simulations 0] and other meth- 
ods §§. 

One way to go beyond this minimal model is to take 
into account the fact that intra- and inter-plaquette ex- 
change integrals might have different values, say Ji and 
J[ , in particular because of the distortion of the lattice 
described in Ref. j2|. However, it seems difficult to ex- 
plain the magnitude of the gap without assuming that 
the ratio J[/ ' J\ is unphysically small. 

In fact, in analogy with the non-depleted square 
lattice 1 10 , it has been suggested that a spin gap 
can be opened by frustration e.g. through the inclu- 
sion of exchange integrals between next nearest neigh- 
bours |U|,||,[ll|], a very reasonable assumption as far 
as quantum chemistry is concerned. This model corre- 
sponds to the Hamiltonian: 



H — J\ S Si.Sj + J2 ^ Si.S 



J 3 1 



(1) 



where J\ (J2) is the exchange integral between (next) 
nearest neighbours on the lattice shown in figure [lj. 
Again, to account for the detailed properties of CaV^g, 
it might prove necessary to allow for different values of 
intra- and inter-plaquette exchange integrals, but these 
differences can be neglected in a first approximation. 
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FIG. 1. Schematic structure of the V4O9 
oxygen have been omitted for clarity. 
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A number of results have already been obtained on 
that model [|6p|J^,[ll[| . According to the "spin- wave" -like 
calculation of Starykh et al. |B[, there is a phase with 
a spin gap and plaquette order when 0.25 < J2/J1 < 
0.8. The cluster expansion of Gelfand et al. || suggests 
that the model has a plaquette-like ground state when 
J2/J1 = 1/2 with a large gap A/Ji w 0.5. This value 
of the gap was confirmed by White ]Tl| ] using a density 
matrix renormalization group (DMRG) calculation, who 
also showed that the gap should open for a frustration ra- 
tio of J2/J1 ~ 0.05, i.e. earlier than predicted by Starykh 
et al. Even if these studies are somehow convergent, a 
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number of points remain to be definitely clarified. First, 
the nature of the ground state in the intermediate frus- 
tration domain has only been determined by perturbative 
methods. In principle, they are valid when the pertur- 
bation involves small parameters, which is not the case 
here. More importantly, the values of the coupling con- 
stants could not be deduced from experiments so far. In 
particular, the interpretation of the temperature depen- 
dence of the susceptibility is not conclusive ||]. Alter- 
native ways of getting information about the exchange 
integrals directly from experimental results are clearly 
needed. 

In order to address these issues, we have performed ex- 
act diagonalizations of small clusters with 8, 16 and 32 
spins with periodic boundary conditions. To reduce nu- 
merical effort, we have taken advantage of all the symme- 
tries of the Hamiltonian: The translations (the basis vec- 
tors are shown in Fig. |l]), the point group C4 (the centre 
of rotation is at the middle of a plaquette) and the spin 
inversion. The elementary cell has four atoms labelled by 
a = 1,2,3,4 (see Fig. [l]) . In the classical limit, the system 
will exhibit Neel order for J2/J1 < 1/2 and a collinear 
order with alternating rows (or columns) of up and down 
spins for J2/J1 > 1/2. A state with Neel (resp. collinear) 
long-range order will then be defined by its wave vector 
(tt,tt) (resp. (0,7r)) and by the orientation e a = ±1 of the 
spins of an elementary cell (ei = —62 = £3 = —£4 = 1 for 
Neel order, t\ = £2 = —£3 = —(-a = 1 for collinear order). 

To test the presence of magnetic long-range order in 
the system, we have calculated the ground state energy 
En, the singlet-triplet gap A at, the magnetic suscepti- 
bility xn = 1/(NAn), and the staggered magnetizations 
Mn(Q) corresponding to Neel and collinear orders and 
defined by: 



m n(Q) 



1 



N(N + 2) 



($| ^e a e^4J |$), (2) 



where t\ is the position of the centre of the plaquette i, 
and e a depends on the nature of the phase. The nor- 
malisation of the staggered magnetizations is chosen so 
that the order parameter is independent of the size in 
a perfect Neel or collinear state. In an ordered phase, 
these various quantities should have the following finite 
size scaling O]: 
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M N (Q) = M*(Q) 
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Actually, for J2 7^ 0, these asymptotic laws are not 
very well satisfied for the available sizes (8, 16 and 32 
sites) p3[ , and the information we could extract on that 
problem is only qualitative. The ground-state energy has 



always a scaling reasonably well described by Eq. (3), so 
it cannot really help deciding whether a state is ordered 
or not. Finite size effects for the triplet gap and for the 
magnetic susceptibility, for such sizes, become large when 
J 2 7^ so that no reasonable scaling could be performed. 
However, the scaling of the staggered magnetization was 
satisfactory enough so that we were able to extract use- 
ful informations. Let us start with Neel order for the 
non-frustrated model with intra- and inter-plaquette ex- 
change integrals J\ and J[ studied by two of us in a 
previous paper |3| . We are now in a position to improve 
the finite-size analysis by considering the results for 32 
spins. It turns out that the scaling of Eq. (4) is not yet 
satisfied for such sizes. Meaningful results can neverthe- 
less be obtained by adding further corrections of order 
1/N. They are shown in figure |[ 
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FIG. 2. Staggered magnetisation of the non-frustrated 
model for N=8, 16,32 and extrapolated values. Quantum 
Monte-Carlo results are from Troyer et al. ||. 

They are in reasonable agreement with the Quantum 
Monte Carlo results of Troyer et al. ||. The slight dis- 
crepancy for small J[ / J\ shows however that such a scal- 
ing should not be taken too seriously at a quantitative 
level. For the model of Eq. (1), the Neel staggered mag- 
netization deduced from a similar finite size scaling anal- 
ysis is shown in figure ||. It vanishes for J2/J1 ~ 0.2. 
This is larger than the value J2/J1 — 0.05 deduced from 
a recent density matrix renormalisation group (DMRG) 
calculation by White |ll[] , which again suggests that large 
clusters have to be studied to get quantitavite estimates. 
As far the collinear order is concerned, we first note that 
the 8 site cluster cannot be used because the collinear 
order is frustrated by periodic boundary conditions in 
that case. Using Eq. (4) with 16 and 32 sites, we found 
that there is a non-vanishing collinear order parameter 
for J2/J1 > 0.7, in reasonable agreement with Starykh 
et al. This order parameter drops abruptly, which 
suggests that the transition with the disordered state is 
first order, as in the case of the non-depleted lattice p0|. 
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So, our results are consistent with the previous results 
that there is no magnetic long-range order for intermedi- 
ate values of J2/J1, although very precise bounds cannot 
be deduced from a finite-size scaling of the results for 8, 
16 and 32 sites. 

Exact diagonalizations turn out to be very useful to 
study the nature of the intermediate phase. The basic 
idea is that, for intermediate values of J2/J1, the sys- 
tem will more or less behave as isolated plaquettes. For 
J2/J1 = 1/2, this can be understood very simply for 
the following reason: Let us write the Hamiltonian as 
H = Hp + Hip where Hp describes independent plaque- 
ttes while Hip describes the coupling between them, and 
let us denote by |$o) the ground state of Hp. Then, for 
J2/J1 = 1/2, (®q\Hip\& ) = 0, which shows that the 
perturbation due to Hip will have a very small effect on 
l*o). 
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FIG. 3. Staggered magnetisation and collinear order pa- 
rameter for the Ji — J2 model for the N=16,32 sites clusters 
and extrapolated values (including also N=8 for the staggered 



In order to check this assumption, we have computed 
(Si-Sj) on different links of the lattice: Along the side of 
a plaquette (PI), along the diagonal of a plaquette (P2), 
or along inter-plaquette links (IP1 and IP2 for and 
J r 2 respectively). A signature of such a phase would be 
a set of correlation functions (Si.Sn) close to their val- 
ues in the pure plaquette phase Jl3. Results for N=32 
and N=16 are shown on figure |4[ The general features 
of these correlation functions are the same for both sizes. 
There is an abrupt variation of all the spin correlations 
which occurs for J2/J1 ~ 0.67 for N=32 and J2/J1 ~ 0.6 
for N=16. This difference in the "critical" frustration is 
the main finite size effect. These curves give us infor- 
mation about the nature of the disordered phase: For 
0.2 < J2/J1 < 0.7, the correlations on all the links are 
close to their values in the plaquette phase. The best 
agreement is for J2/J1 = 1/2. For this particular value 
of the frustration, the correlation functions on the dif- 



ferent links are comparable to their values in the pure 
plaquette state: -.4623 vs -0.5 for PI, .2201 vs 0.25 for 
P2, -.0963 vs for IP1, and -.0009 vs for IP2. Hence, it 
seems clear that the ground-state is a plaquette resonat- 
ing valence bond state. Note that there is no sign of the 
ground-state being degenerate, which allows one to elim- 
inate the alternative broken symmetry state proposed by 
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FIG. 4. (Si.Sj) on different links for the N=16 (a) and 
N=32 (b) clusters: PI (•), P2 (■), IP1 (o) and IP2 (□). 

We now discuss the low-lying excitations in the inter- 
mediate frustration regime where the system has a pla- 
quette ground state. For a single plaquette of four spins, 
the ground state is a singlet of energy — 2Ji + J2/2 as long 
as J2 < J\. However, the first excited state is a triplet 
only if J2 < Ji/2. Beyond that value, there is a singlet 
state in the singlet-triplet gap. Let us see what remains 
of that picture when plaquettes are coupled. The lowest 
energy levels versus J 2 / J\ for the N = 32 spin cluster are 
shown in figure ||. The other clusters have the same be- 
haviour. The ground state is always a singlet totally sym- 
metric under the symmetries of the Hamiltonian. The 
nature of the first triplet excitation depends on the value 
of J2I J\- When J2 / J\ is small, the first triplet excitation 
has a wave vector Q = (ir, ir), while for J2/J1 > 0.56 the 
first triplet has a wave vector Q = (0, ir). This change in 
the position of the lowest triplet in the Brillouin zone is 
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consistent with the prediction of Starykh et al. || and of 
Gelfand et al. || who both found such a change around 
J2/J1 = 1/2. A Schwinger boson mean field calcula- 
tion |n| actually suggests that the position of the maxi- 
mum of the spin structure factor moves continously from 
(■7T, 7r) to (0,7r). Now, more importantly, we found that, 
for intermediate values of J2/J1, there are singlet states 
in the triplet gap when 0.49 < J2/J1 < 0.7. One of these 
excitations has the same wave vector as the ground state 
but has a d-wave symmetry. For the N — 16 clusters, 
this energy level crosses the ground state energy but this 
is probably an artifact of this cluster. The existence of a 
first excitation which is a singlet had not been reported 
so far for this model. 



been performed leading to J 2 ~ 530K. The local geome- 
try being essentially the same in these compounds and in 
CaV40g, the exchange integrals should be comparable. 
So these results suggest that J 2 should be of order 500 K, 
and that J2/J1 should be around 0.7. Note that large val- 
ues of Ji and J2 are compatible with the reported value 
for the gap (107-ftT) because the gap is much smaller than 
Ji close to the boundary to collinear order. Observing 
singlet states in the singlet-triplet gap would be a direct 
confirmation of this picture. 

We acknowledge very useful discussions with L. Levy 
and C. Lhuillier. The numerical calculations reported 
here were possible thanks to computing time made avail- 
able by IDRIS, Orsay (France). 
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FIG. 5. Lowest energy levels for N=32. Solid symbols are 
for spin state, while open are for spin 1. Dashed lines in- 
dicate the position of the change of the nature of the first 
excitations described in the text. 

Testing experimentally the presence of singlet states 
in the singlet-triplet gap in CaV-iOg should give useful 
information about the frustration ratio for the follow- 
ing reasons. We first note that two related compounds 
have been studied recently, suggesting that the value of 
the coupling constants could be quite different from the 
commonly accepted value, i.e. J 2 j ' J\ ~ 1/2 and J\ rang- 
ing from 100K to 200K. A parent compound of CaV^g, 
CaV307, which has a similar structure as far as the local 
geometry in the VO plane is concerned, has been studied 
by Harashina et al. [^]. Using neutron diffraction, they 
have shown that it has a collinear magnetic long range 
order. This is the signature of a large frustration: A 
modified spin wave calculation of Kontani et al. Jl8| sug- 
gests that this state is the ground state only if J 2 /Ji is 
bigger than 0.7. More recently, a quasi-one dimensional 
vanadate, NaV20s, has been studied p9[ . This mate- 
rial consists of spin 1/2 chains corresponding to corner 
sharing VO5 square pyramids. So it involves only one 
coupling constant which should correspond to J 2 . A fit 
of the temperature dependence of the susceptibility has 
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